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ABSTRACT 


The interaction lengths of pions and protons in iron have 
been measured using an ionisation spectrometer composed of 
alternating layers of iron and plastic scintillator. These 
measurements cover ;in energy range from 9.3 to 18 GeV. The 
interaction lengths ware determined by accurate statistical 
analyses of the experimental data using the maximum-likelihood 
method. The dependence of the interaction length on the para- 
meters used to define an interaction was studied, and the results 
reported employ parameters chosen to minimize the percentage 
uncertainty in t’-e Jntciactlon length. The mean interaction length 
of pions was found to be app cox - * lately 2 0% greater than that of 
protons. No significant dependence on the energy of the incident 
primary was found in the range measured. Consistency was found with 
the great variety of accelerator measurements made to date, all of 
which report greater interaction lengths than those reported by the 
cosmic ray measurements . 



INTRODUCTION 


As energetic charged particles traverse natter they lose 
energy in small amounts via an almost continuous process of 
excitation and Ionization of the atoms and molecules of the 
aboorber material. However, strongly interacting particles 
(hadrons) also have n probability of interacting with in situ 
nuclei. The strong interaction is characterised by the produc- 
tion of one or more additional strongly interacting particles 
(usually tt or K mesons) Mid hence constitutes a much larger 
and discrete energy loss by the incident hadron than is typical 
of the more frequent atomic collisions. 

The interaction probability of hadrons in hydrogen has been 
studied extensively, both as a function of the. incident energy 
and the type of incident hadron. The probability for interactions 
in materials other than hydrogen has been leer, extensively studied. 

In this paper we report the results of measurements on the. probability 
for strong interactions of incident pions and protons in iron. 

Similar measurements have been reported previously, but these 
earlier works suffered from large uncertainties in the incident 
energy or from less accurate methods of identifying the strong inter- 
action. In this experiment, lnomen turn- analyzed beams of pions and 
pro tons produced by the Alternating Gradient Sychrotron were employed. 
The de Lector was composed of modules consisting of alternating layers 
of lion and plastic scintillator. The output signal from each module 
was digitized and recorded on magnetic tape, so that each event could 
be reconstructed end analyzed In terras of a variety of parameters. 
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In the present experiment, ae well as in most of the 
previous work, the parameter that is measured is the module 
number of the detector in which the interaction took place. 

If the probability of an interaction is independent of the 
particle energy, then the probability of an interaction per 
unit thickness of absorber is independent of the depth in 
the material. Hence the probability of a particle traversing 
a depth- £ without interacting is given by 

POO - c" £/x . (l) 

The constant X is the mean interaction length, and is related 
to the total interaction cross section, o, by 



crA v p 


where A and p are the. average a toraic muss and the density of the 
target material, and A is Avagadro's number. The parameter X 
is determined from the distribution of interactions in the detector. 

Details of our experimental technique, procedures for data 
reduction and analysis, and comparisons with previously reported 
results. are described in the following sections. 

APPARATUS 

The cosmic ray detector used in this experiment is designed to 
measure the energy spectra, charge composition, and arrival direc- 
tions of cosmic rays in the energy range from 10'*"® to lO^ 1 eV« The 
data for this experiment were accumulated during a calibration run 
performed in. the G10, +4.7° test beam at the Alternating Gradient 
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Sychrotron at Brookhaven National Laboratory. The instrument, 
show schematically in Fig. 1, contains a particle identification 
section which serves to define the charge and trajectory of the 
incident particle and an ionization spectrometer to determine 
the energy of the incident particle. 

The spectrometer consists of 12 tungsten modules and 7 iron 
modules. For a portion of the present measurements , 4 additional 
iron modules were incorporated in the spectrometer, using elec- 
tronics irora 4 tungsten modules. 7. he spectrometer is designed to 
permit unambiguous identification of incident electrons through 
the rapid development of electromagnetic cascades in the tungsten 
portion of the. detector. Each tungsten module, consists of a 
tungsten plate 6.3 g /cm thick (0-37. cm or approximately 1 radiation 

n 

length) and a scintillating plastic sheet 0.6b g/ci:r thick (0.64 cm). 
Each of these modules is viewed separately by a pair of photomultiplier 
tubes . 

The spectrometer is also designed to measure the energy' of primary 
cosmic ray protons principally' through the response of the iron 
modules to energy' deposited by nue 1 csr-c lec tromagne tic cascades. In 
the detector used in the present experiment, the arrangement of 
absorber material (iron) and active elements (plastic scintillator) 
is such as to optimize the energy resolution of the spectrometer by 
minimizing fluctuations in its energy response. Conventional ioni- 
zation spectrometers used in cosmic ray experiments exhibit an energy 
resolution which is dominated by two major sources of fluctuations. 
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One is caused by uncertainty in the location of the first 
hadronic interaction, which leads directly to an uncertainty 
in what fraction of the primary energy escapes out the bottom 
of a finite - depth detector. The other fluctuation results 
from infrequent sampling of the. energy deposited in electro- 
magnetic cascades which are initiated by the neutral pi ons produced 
in high energy hadronic interactions. In the present experiment, 
each iron module consists of 66.4 +0.7 p/cm" (approximately one- 
half of a nuclear interaction length) of iron and is viewed by a 
pair of photomultiplier tubes. This enables location of hadronic 
interactions to within one half a r .clear interaction length. The 
iren in each module is arranged in layers interspersed with three 
0.64-ca thick plastic scintillator sheets. The active elements are 
thus located approximately every 1.5 radiation length of absorber 
material so as to thoroughly sample electromagnetic, cascades. This 
spectrcraerer design has greatly improved the energy resolution over 
previously reported results."*® 

The output signal of each of the tungsten and iron modules was 
calibrated using cosmic-ray muons, and non-interacting protons 
produced by the accelerator. Calibration runs were taken both before 
and after the proton and pio-.r data were accumulated. The response of 
each of the modules to equivalent single relativistic muons was 
determined by pulse height analyzing events in which energetic 



- 5 - 


cosrnic ray muons were incident on the detector. Such events 
were selected by demanding that alternate modules situated before 
and after each of the modules being analyzed had pulses corres- 
ponding to the passage of a single relativistic particle. 


ANALYSIS OF DATA 

In the approximately four nuclear interactions lengths con- 
tained in the ionization spectrometer (six Interaction lengths 
after the additional modules were added) almost all incident hadrons 
participated in a nuclear interaction. The problem was to determine 
reliably in which module the interaction occuired. It was assumed 
that at least one additional high energy particle (usually a pica) 
is; created in any nuclear in! erection. > 1 1 u n the mod 'Jos located at 
greater depth than the in to r a cl ion p.hnw an increase in signal associated 
with the passage of more tlum one* relativ is t.ic particle. When a 
neutral pion id created in the nuclear interaction, the resulting decay 
to two gamma rays and subsequent electromagnetic cascade produces rainy 
ionizing particles in the following modules. 

In order to determine the module in which the nuclear interaction 
took place, the* pulse height f ro,« each module is divided by the average 
equivalent muon pulse height. The numbers thus obtained, which are 
assumed to represent the number of icrizing particles in each module, are 
used to determine the interaction location. 

To illusfiate the nature of the problems encountered, plots of 
the number of ionizing particles as a function ol. module number are 
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ahown in Fig. 2a and 2b. Figure ?.a illustrates a case, which 
ia easy to analyze. A single incident particle passed through all 
12 tungsten modules end two iron modules. The outputs of these 
modules show only the usual statistical variations associated with 
the passage of one relativistic particle. In Iron Module 3 (Fc 3) 
an interaction occurred which produced one additional particle. 

The subsequent nuclear and elect .vomaguetic cascade is seen to build 
up in the deeper modules. Figure 2b shows a more typical cane arsd 
is clearly more difficult to analyze unambiguously. In this case 
two isolated nodules (Tungsten 1 end 11) recorded pulse heights consistent 
with the passage of several relativistic ionizing particles but are 
followed by modules with outputs corresponding to one relativistic 
particle.-. It is av.'-u -'eci thr.t the .so isolated high-output signals 
are due to stocaotic variations in the nunher of scintillation photons 
detected, photomultiplier tube noise, or the production of delta rays 
that stop in the detector. This conclusion is supported by the 
observation that these "spuriously" high outputs appear with the same 
frequency even when the. incident particles are muons. In the Fe modules 
the effect of Landau fluctuations is greatly reduced by the averaging 
of the pulse height from the three scintillators. In nil of the 
analyses reported here, at least two consecutive modules indicating 
more than one relativistic ionizing particle were required to define 
an interaction. This effectively requires that charged pions be 

n 

produced with a range of wore than 66 g/cm^ in iron or that a tr° be 
produced with enough energy (~ 600 KeV) ro that the associated electromagnetic 
cascade penetrates more than 9 radiation lengths of iron. 



y; 
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In order to analyze the 5 x 10^ events recorded during the 
course of this experiment, a computer program was written to 
process the data. The data recorded from each incident particle 
were searched for two or more consecutive modules with outputs 
above a threshold level. The interaction was defined to have 
started in the first of the modules in which the threshold was 
exceeded. Hot.li the value of the threshold and the number of consec- 
utive modules required to exceed this value were adjustable para- 
meters. After reduction of the data for a particular run, a summary 
of the information for each module was printed. This summary included 
the total number of events analyzed, the number of times a single 
particle, passed through each module, rl.o number of tines the signal 
was spurious: (i.e, the signal exceeded the threshold value but the 
signal in the next nodule did notv , and the number of times n particle 
interacted in each module. 

In st 'tie instances spurious effects generate a high signal in the. 
modul: just before the one in which an interaction occurred. Tsis 
spurious signal causes the wrong interaction module to be selected. 

If all modules were spuriously high, the same percentage of the time, the 
slope of the interaction curve would not be affected. However, this was 
not the case; the spurious rate varied from much less than 1% to slightly 
more than 5%. The number of interactions found in each plate was 
therefore corrected, as follows, to account for the affects of spuriously 
high events. 

If module i produced a spuriously high signal a fraction of the 

i 

It 
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tiiae S^, then the observed number of events in module i, should 
be reduced by S. H,. because of the spurious signals. However, the 
number of events in module i. should be increased by toe mmber of 
events falsely assigned to the i-1 module. Thur the correct number 
of events observed in module i is given by 

Hi » Hj. (1-Si) + S i-1 ( 3) 

The mean interaction length w«i» found by ueueiruinlug the slope 
of the exponential decay of the number of interactions in each of 
the iron modules. In all of these determinations the interactions 
in both the first and last iron modules were ignored. The first iron 
module *• ***> ignored because of predicted bach sc.a ttering effects,^ and 
mote imp ' ly, because this module is not symmetrically embedded in 

iron. The showers created by interaction in the dteper portions of 
this module are detected by only one or two scintillator sections, and 
are hence seen less efficiently than interactions occurring in the front 
of the module. Such low efficiency everts will often be registered 
as having started in the subsequent Iron module, but there is no corres- 
ponding nodule in front of the first one to add to events detected etr. 
starting in the first, module. Thus the efficiency of the first module 
is predictably too low. The interactions detected in tire last module 
were also ignored because of a lack of symmetry in the bad-scattering 
and because there was no effective way to determine the number of times 
the last detector was spuriously high. 

STATISTICAL ANALYSIS OF THE EXP ®|r MENTAL DATA 
The processes by which energetic hadrons interact in an ionization 
spectrometer are stochastic, and the relationships describing the 
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expected response of the detectors are statistical. In the present 
work, th- method of maximum likelihood was used to determine inter- 
action lengths from distributions of interactions in an ionization 
spectrometer. This ne th orl is readily applicable because each 
pai title has a well defined a. priori^ probability of interacting as 
a function of detector module number, independently of the inter- 
action probabilities of each of the other particles. 

Another approach, sometimes used to determine the interaction 
length, < .'--.ploys a less cumbersome technique. The logarithm of the 
observed in Lei action- d f s Ir i bution is fit by a weighted least-squares 
to a straight line. The applicability of the least-squares method 
relics on the indope’idt-nco of the numbers of particles observed t 
interact in each deter! or nodule and the appropriate choice of tub weighting 
function. In any attenuation experiment such as reported here, the 
numbers ci particles interacting in the different detector mode Les 
are not independent, and are in fact related by the expression, 
k 

N = £ n , ' (A) 

where N is the total number of interacting particles, k is the number 
detector modules, and the n^ are the numbers of particles observed 
to interact in each of the raouules. That is, the probability of interacting 
nodule i is coupled with the probability of not interacting in the 
i-1 proceeding modules. In the present experiment, the total detector 
thickness is large, approximately 4 interaction lengths of absorbing 
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material. Therefore., the total number of interacting particles 
is a large fraction of the total flux through the detector, 
fluctuations in the number of particles interacting in any module 
propagate as fluctuations in the flux through succeeding detector 
modules. Thus, the approximation that the are independent cannot 
be justified. Its use would lead to results which are influenced 
by coupled statistical fluctuations in the numbers of interacting 
particles in successive detector modules. 

For each data set in the present experiment, an interaction length, 
X, was determined by maximizing the likelihood, as a function of X, 
given by the product of all the a pri ori probabiiivi.es for each event 
in the data set, 

T1 = probability that a particle, present at.J!=0, will arrive at 
- o / \ 

t is e 5 where X is the nuclear interaction length. The probability 

that a particle at £ will interact between £ and £-Hl£ is dX, Since 

X 

these probabilities are independent, the probability that a particle 
present at £~-0 will arrive at £ and will interact between £ and £+A£ 
is given by 

P(£,A£,X) ~ f e d/X, or (5) 

P(£ A£ X) = e~ l/X (1 - e~ hZ/X ), (6) 

where A£ is a module thickness, and £ is total thickness of modules 
up to the one in which the particle interacts. However, in this experi- 
ment, the only particles selected are those which interact between 
£=0 and £=L, where L is the total thickness of all the modules. Thus , 
the probability of interacting in the whole detector, treating it as 
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one module, Is (from Eq. 6) 

P(0, L, A) => 1 - e" L/x . O) 

So the probability of interacting between SL and &+A2. for all 
particles which arrive at £-0 and which do interact between 
SL- 0 and £=>L is given by 

P(£ : A?,, L, X)“ C.~ l/X (1 - e" A£ - /X ) / (1 - e -L/?1 ). (8) 

An analogous expression is commonly used in bubble chamber experi- 

1 ^ 

nier.ts to determine particle lifetimes. 

For large L/A and small Af/A, Eq. (8) reduces, to the well known 
limiting expression 


P(?,, At, L, A) 


M. _-£/X 


lim L/A-> « 
A&/A-+- 0 


(9) 


However, since the limit criteria are not satisfied in the present 
work, Eq. (8) is the appropriate expression for the a p riori proba- 
bility end is used to calculate the likelihood function. 

The likelihood function is given by 


TT P j _ (£, At, L, A). 
i“l 


( 10 ) 


Using P^ from Eq. (8) gives 

X - |] (1 - n~ Ll/X ) / (l-e~ L/A ), (11) 

i=l 

where i indexes each of the N independent events in the data set. Tire 
interaction length, to be determined from the data set, is that value of 
the quantity. A, for which the likelihood, , fa a maximum. Since ^ is 


- 12 - 


a positive definite quantity, it is a maximum when 

***-> 

- Z |l£i+ ln(l-e~ 6i/X ) -in(l-e“ L/X )| (12) 

i=l L x J 

is a maximum. Analytically, in(£) C3n be maximized by setting its 

first derivative with respect to X equal to zero. The resulting 

transcendental expression for X must be solved by numerical techniques. 

In the present work values of ?.n (£) were computed for each data set 

for a range cf values of X. An ( £) was found to be nearly normally 

distributed in reciprocal X, 1/X, for all data sets, independent of 

size. For data sets containing more than 2500 events, An(£) was found 

to be normally distributed in X as well. The observed behavior of the 

likelihood function in X and reciprocal X has been predicted for this 

type of data in Review of Particle Properties (1970) by Barbaro-Gal tieri ' 

3 3 

et al., ' and can be derived from the analytical behavior of Kq. (12). 

The measured uncertainties wore determined from the numerical 
calculations from that change in X, 6X , for which P.n( t £) decreases 
by the amount 0.5 from its maximum value. This corresponds to the 
one standard deviation points of a normal distribution. A limit on the 
best accuracy obtainable in X fora total number of events, N, for a smalx 
AA/X and for. a finite L/X has been presented in Ref. (13). In the present 
work, in which both ASt/X and L/X are finite, the minimum uncertainty 
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in X can be shown to be 


6X 


min 


X f e AS./x / M-/X \ 2 

7¥l | e A *A-l / 



( 13 ) 


All standard deviations determined in the present work slightly exceed 
<5 X „ -■ _ . The values of X and 6X determined from the analyses .just 
described are presented in the next section. 


kESULTS 


Because the analysis program was written in a general manner, 
it was relatively easy to determine both the slope of the interaction 
curve and the quality of the fit as a function Of numerous parameters. 
The data were analyzed as a function of the number of consecutive 
"high" modules required to define an interaction, as a function of the 
number of equivalent muons needed to define an interaction, and as a 
function of different numbers of muon equivalent particles in the 
first and second interaction modules. 

The slope of the interaction curve was calculated for several 
different runs using the requirement of outputs exceeding a threshold 
from two or throe consecutive modules to define an interaction. 

For eacli run the slopes determined using the requirement of two 
consecutive modules agreed well within the experimental uncertainties 
with tlie slopes determined using the requirement of three consecutive 
modules. However, the statistical accuracy associated with the three 
module measurements was somewhat less. The data from one less module 
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could be used since the number of interactions occurring In 

the next-to-last module obviously could not be determined with 

the same requirements placed on the other modules. The loss 

of information which might have been provided by this one 

additional module increases the statistical uncertainty. Since 

the determination of the slope was independent of whether the 

two or the three consecutive-module, requirement was used,, and because 

greater accuracy was achieved with the two module requirement, all 

succeeding determinations v arc nn.de with the requirement that two 

consecutive modules had outputs which exceeded the selected thresholds. 

The threshold chosen to define an interaction selects the 
physical processes to be observed . As the threshold for determining 
an interaction is raised, interactions which have low multiplicity 
and also fail to develop a significant shower are not detected. 

However, interaction; with high multiplicity will continue to be 
detected. Most interactions in which one or more rr°'s are produced 
will also be detected due to the rapid build-up of the electromagnetic 
shower following the decay of the rrO into two gamma rays. At the 
energies investigated in this work, higher multiplicity events have 
lower probabilities. Since increasing the threshold requirement will 
bias against some low multiplicity, higher probability events while 
continuing to select the high multiplicity, lower probability events, 
the interaction length is expected to rise as a function of increasing 
threshold requirement. Figure 3 shows a plot of the measured inter- 
action length for both 9.3-GaV protons and positive pions as a function 
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of the threshold. The interaction length is seen to increase 
by about 20% as the threshold requirement is raised from 2 to 
7 single particle equivalents. The measured interaction length 
seems to remain roughly constant or decrease slightly for still 
higher thresholds. However the statistical accuracy for the higher 
threshold requirements rapidly becomes worse because not many events 
satisfy the lather unusual threshold requirements. Tor very low 
threshold requirements the. statistical uncertainty in the determina- 
tion of the slope increases for similar reasons. For example, as 
the threshold is lowered to one ionizing particle the interaction will 
almost alvaye be detected in the first few tungsten modules and no 
useful data for the determination of the slope in the iron moduelc 
wi 1 1 be avu i J ah le . 

Not much data arc available for the particle multiplicities to 
be expected from nuclear interactions in iron at the ■ .’.ergies inves- 
tigated in the present wor-k . Thus the interaction mean free path 
was determined for a range, of values of the thresholds. Figure A 
shows a plot of the statistical uncertainty in the meon free path 
as a function of threshold. For these data there is a very clear minimum 
in the percentage uncertainty for a threshold requirement of 3.0 muon 
equivalent particles. The other data (13.8 and 17.8-OeV protons ar.d 
9.3-GeV pious) also show a minimum in the uncertainty at a threshold 
value of 3.0 inuon equivalent particles although in these latter cases 
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cases the minimum is not so pronounced. This minimum in the 
uncertainty occurs even though tiie total number of interactions 
found in the iron modules is a monatonically increasing function 
of the threshold value until the threshold is approximately 8 
muon equivalent particles. The uncertainty in the determination 
of the slope depends both upon the statistical accuracy of the data 
and upon how well the data fit the assumed exponential decay function. 

The foot that for values of the threshold greater than 3 the quality 
of the fit becomes worse while the total number of interactions 
detected increases, indicates that data for those higher thresholds 
do not fit the hypothesis of a simple exponential decay so well. The 
threshold which resulted in the minimum statistical uncertainty in 
X was chosen. 

Once criteria for the number o£ modules (2) and the threshold value 
(3 muon equivalent particles) had been determined, the effect of adjusting 
the threshold of the second module while holding the threshold in the 
first at 3 was investigated. For some cases, particularly for inter- 
actions which initiate electromagnetic cascades, the signal level in the 
second module is higher than in the first. It was thought that a higher 
threshold for detection in the second module would more nearly represent 
true events and be a better criterion for event selection. Tire data for 
9.3-GeV protons were analyzed by varying the threshold for the second 
module from 1 to 5 muon equivalent particles while the threshold for 
the first module was held constant at 3 muon equivalent particles. As 
can be seen in Fig. 5, the minimum uncertainty in the determination of 
the slope occurs with the same threshold for both the first and second 
modules. While adjusting the threshold value in the second module changed 
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the accuracy of the dope determination; no appreciable change 
In the value of the Interaction length was detected. Thus all 
final results on the determination of Interaction lengths presented 
here are based on the requirement of two successive modules with a 
signal output above 3 times the single particle level to define an 
interaction . 

The results of the determination of mean interaction lengths 
for protons and pions arc presented in Tables I and II. The number 
of modules employed, the total number of interactions analyzed, 
and the measured mean nuclear interaction length are given for each 
incident energy. The uncertainty in the interaction length determined 
by the maximum likelihood method is also given, along with the theore- 
tical minimum statistical uncertainty based on the total member of 
interactions detected. Fach of the experimentally determined interactions 
lengths coniains an additional uncertainty of +1% clue to the uncertainly 
in the mean density of the detector material. 

It needs to be strongly stressed that in Tables I and II the values 
of the interaction lengths and the. associated uncertainties are given 
for our fcest-fit definition of an interaction (3 muon equivalent 
particles in each of two consecutive modules) . As Fas been discussed 
previously, the value of the interaction length can be changed con ,derably 
outside the stated errors by appropriate changes in the definition of an 
interaction. 

The mean Interaction lengths are observed to be independent 
of the number of iron modules employed to within the statistical uncer- 
tainties in the measurement. No significant dependence of the interaction 
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length on energy is detected. On the hypothesis that the inter- 
action length is independent of energy V = 0.09 for protons and 
P " 0.95 for pions. The mean interaction length for pions is 
approximately 20% greater than that for protons. 

Some measure of the overall accuracy of the determination of 
the mean interaction length can be obtained from a comparison of the 
experimentally determined uncertainties with the theoretical minimum 
uncertainties predicted from Eq . (13). The fact that, the experimental 
uncertainties are only slightly larger than the theoretical minimum 
indicates that the assumption of a single exponential decay function 
is reasonable. Departures of the data from the assumed exponential 
decay could be caused by variation of the interaction length with 
energy, unequal efficiencies of the detector modules, and contamination 
in the be a"). 

The effect of energy dependence on the determination of the mean 
interaction length, is expected to be very small. The incident particles 
lose less than 10% of their energy in traversing the entire detector, 
even for the lowest energy (9.3 GeV) protons. Since the interaction 
probability is at worst a very slowly varying function of the energy, 
the small change in incident energy should not appreciably affect the 
determination of the near, interaction length. 

Contamination of the proton beam by other hadrons was small. 

v 

Efficiency for rejection of ir and K mesons by the gas Cerenkov detector 
was measured to be better than 98%. Since the flux of protons always 
represented more than half of the total flux, the contamination in the 
proton beam was always less than 2%. 
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There was, however, some K meson contamination in the tt 
meson beam. At both 5.3 and 17.8 GeV the ratio of IT*" to ir + Q,; 
the production target was 12%. Because of the long flight path 

t l 

(57.35 m) and the lower energy of the K , about half of the K 

4 * 

mesons decayed before reaching the detector. The ratio of K 
to 'it + incident upon the detector was 5.1 and 7.8% at 9.3 and 
17.8 GeV respective] y . The correction for the K 1 " contamination 
in the beam has been made assuming the ratio of the cross sections 
0(K + ,F > to o(iT + I ,Fe) will be the same as the ratio of the cross 

/ -h 

section o(K' , D) to a(r , 11) of approximately 0.71 given by Denisov 
15 

et al. The correction for the It meson contamination has been included 
in the results presented in Tabic II. The corrected value for the. 
t: + interaction length in iron at 17. C CoV/c. of 164 +9 g/cai^ is in 
excellent agreement with the value of the r interaction length measured 
at 20 GeV/c of 162 -VI g/cir. obtained by interpolating between the values 
given for Al and Cu by Allabv et al.' L ® Karrie et al.^ have reported 
a mean Interaction length of 130’ g/cm?- for 3.75 GeV/c iri in iron. This 
determination is based on the interaction criterion of scattering 
of the incident pion, so thaC the much higher interaction cross section 
may be due to an effectively lower threshold criterion. 

The effect of detector efficiencies upon the determination of the 
mean interaction length has been investigated, first, as a check 
for variations, in the efficiencies for detection of an interaction, 
the total number of interactions detected in each module for all proton 



runs was plotted as a function of module number. For each module 
included in the data analysis the sut; of the interactions differed by 
no more than ?.% from the expected fit to the data. We thus conclude 
that all detecLors have the same efficiency for the detection of an 
interaction to within a few percent. Secondly, the muon eqtiivalent 
particle calibration for the first useful iron module was arbitrarily 
changed by 10/,, a number several times the assigned uncertainty. The 
value of the mean Interaction length changed by less Ib.m 1%, but the 
quality of the fit became worse. 

Cone lesi ons 

The nean interaction lengths for 9,3-, 13. 8-, and 17 . 8-GeV protons 
and 9.3- and 17. 8-GeV positive pions in iron have been determined. The 
value obtained depends somewhat upon the signature required to define 
an interaction. Exact comparison with the other determinations of the 
mean interaction length is not possible because other experiments have 
slightly different criteria to define an interaction. The results of 
this experiment and those of previous workers are shown in Table III. 
Because of the accurately determined incident energy and large number of 
events, the precision of the present results is generally greater. 

Ko statistically significant variation of the mean interaction 
length for protons or pions as a function of energy was observed in 
this work. However, with only two exceptions, the data presented in 
Table III show a systematic. 5% difference between measurements of the mean 
interaction length made with cosmic rays and those made with accelerator 
produced protons. A correction for pion contamination in the cosmic 
ray flux would lower the value obtained in uncorrected cosmic ray experi- 
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ments by about 5 % thereby Increasing the discrepency. 

All of the rscent accelerator experiments have employed both 
a similar detector constructed of alternating layers of iron and 
scintillator and similar criteria for defining an interaction. Thus 
the agreement among these experiments is not unexpected. The analysis 
of the data in the present work has shown that the value of the mean 
interaction length is dependent upon the criteria used to def.tnp an 
interaction. However most of the changes in criteria from the ones 
finally selected for our analysis tend to raise the value of the mean 
interaction length. None of the various interaction criteria that were 
tried gave a value as small as with those reported by the cosmic rays included 
experiments. The difference may be due to the effectively higher 
threshold criteria in the cosmic ray experiments resulting from less 
frequent sampling and greater module thickness in their detectors. 

Tie determinations of the interaction length using cosmic rays 
Included only particles with energies above 3Q GeV. The observed differ- 
ences in tiie mean interaction length could be explained by an energy 

dependence the interaction cross section. Such an energy dependence 
n 9 

has bven sought by Jones et al. over the energy range 70-800 GeV but 
no statistically significant effect was observed. The accuracy of these 
investigations was, however , not sufficient to rule out an energy 
dependence as the cause of the observed discrepency. 

It is always difficult to make comparisons between experimental 
me a- axemen ts in which different techniques, that may affect the results, 
have been employed. It will be interesting to repeat the mean inter- 
action length measurements using an iron-scintillator spectrometer with 
cosmic rays as a source of protons. 
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Figu res 


Fig. 1 A schematic view of the. high-energy cosmi( — ray 
ionization calorimeter used in this experi men t . 

Fig. 2a Plot of pulse height as a function of module 
number. In this case the. module in which the 
interaction occurred is clearly Fe 3. 

2b Plot of pulse height as a function of 

module number. This plot shows the difficulty 
of making a unique selection of the ini erection 
module . 

Fig. 3 Plot of the measured interaction length as a 
function of threshold requirements. 

Fig. 4 Plot of statistical uncertainty in the determination 
of the mean interaction length as a function of 
threshold requirement. 

Fig. 5 Plot of statistical uncertainty in the determination 
of the mean interaction length as a function of 
threshold requirement in the second interaction module 
with the first held fixed at three equivalent particles. 


Tab Is I 


Mean Interaction Length 


of Protons. in Iron 


Incident 

Number of . 

Total Kunber 

. i 

AA/A in 

Percent 

Energy 

Iren 

of Defected 

. 



(GeV) 

Modules | 

1 

Interactions 

X in g/cirr 

Measured 

Minitaisn 

9.3 

5 

2942 

150.8 +4 ’ 8 
-4.7 

Hr 3 . 2 

+3.0 

13.8 

5 

2804 

151.4 +5 .’\ 
-4.7 

+3,2 

+3.0 

17.8 

5 

3137 

139.0 +4.1 

+2.9 

+2.8 

17.8 

9 

2683 

140.2 * 3 ' 2 
-3.0 

+2.3 

+2.2 
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Table II 

Mean Interaction Lengths of Picns in Iron 


Incident' Number of 


rot 


o 

X in g / c trr 


h\f\ in Percent 


Energy j Iron of Detected i Unccr re rtcd Corrected 


(GsV)' 

Modules 

Interactions j 

for K‘‘ Cents 

nination 

Measured 

Minimum 

9.3 

5 

1291 

166.8 +8,5 

164.5 +3.5 

±5 • 1 

±4 • 4 

17.8 

9 

550 

167.4 +9.1 

163.8 +9.0 

±5.3 

±4.9 





Table III 

Comparison of Experimental Determinations of Kean Interaction 
Lengths of Protons in Iron 


Experimenters 

Incident 
Energy 
(Go V) 

Source of 
Protons 

A in 
g./cr 2 

Comments 

Brenner and Williams ^ 
(1957) 

50 

Cosmic Rayu 

152+7 

Integrated all ha J *-ons 
in cosmic ray flux'- 6 ) 

Ashmore et all c ) 
(I960) 

24.2 

Accelerator 

143+8 

Value interpolated from 
measurements on Al and Cu 

Dobrotin et 
(1965) 

20-600 

Cosmic Rays 

125+16 

Integrated all hadrons in 
cosmic ray flux^ b ) 

Belle ttin.l et all e ) 
(1966) 

19.2 

Accelerator 

119+3 

Value interpolated from 
measurements on C, Al, Cu & Pb 

Andrcnikashvili^ 
: (1968) 

50-1000 

Cosmic Rays 

130+6 

Integrated all hadrons in 
cosmic ray flux^ b ) 

(Vt 

Scnmidt Xo/ 

(1968) 

28 

Accelerator 

143+3 

Minimum statistical error 
+3.3 g/c.m^ 

i \h V. Jones et al.^) 

i U970) ^ 

! Bashindjhagyan at al. ' 

1 (1971) 

I 

20.5-28 

400 

Accelerator 
Cosmic Rays 

139 

132+5 

Combined 20.5 and 28 GeV data 

Corrected for pions in 
cosmic ray flux 

| Kurz et 

! to be published 

i 

18 

Accelerator 

135.3+3.4 

Experimental equipment and 
data analysis similar to 
present experiment 




Table III 



Incident 





Energy 

Sources c-f 

A in 


Experimenters 

(GeV) 

Protons 

g/cn* 

Comments 

L. W. Jones et al.^) 
<1972) 

70-800 

Cosmic Rays 

132.4+2 

Corrected for pions in 
cosmic ray f lvx 

Present Experiment 

17.8 

Accelerator 

139.9+2.5 

Bcnand 3 particles in 2 


successive modules 


(a) 

(b) 

(c) 
'(d) 

(e) 

(f) 

(g) 
00 

(i) 

(j) 
00 


Ref. 1 

Correcting for pion contamination in the cosaic ray flux vcuid lover the value of the mean 
interaction length bv about 52. 

Ref. 2 
Ref. 3 
Ref. 4 
Ref. 5 
Ref. 6 
Ref. 7 
Ref. 3 

Private communication - We are indebted to Dr. Kurz for permission to use these data prior to publication. 
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